Langevin dynamics of the Coulomb frustrated ferromagnet: a mode-coupling analysis. 
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We study the Langevin dynamics of the soft-spin, continuum version of the Coulomb frustrated 
Ising ferromagnet. By using the dynamical mode-coupling approximation, supplemented by reason- 
able approximations for describing the equilibrium static correlation function, and the somewhat 
improved dynamical self-consistent screening approximation, we find that the system displays a 
transition from an ergodic to a non-ergodic behavior. This transition is similar to that obtained in 
the idealized mode-coupling theory of glassforming liquids and in the mean-field generalized spin 
glasses with one-step replica symmetry breaking. The significance of this result and the relation to 
the appearance of a complex free-energy landscape are also discussed. 
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I. INTRODUCTION 

Models with a competition between a short-range or- 
dering interaction and a long-range frustrating interac- 
tion have been recently introduced to explain the slow- 
ing down of relaxation in supercooled liquids and the 
resulting glass transition Q. The underlying picture is 
as follows: any given liquid possesses a locally preferred 
structure which is different than that of the actual crys- 
talline phase, and this local arrangement of the molecules 
in the liquid cannot propagate at long distances to tile 
the whole space and form an "ideal crystal" because of 
ubiquitous frustration. It has been argued]^ that, for 
weak enough frustration, this phenomenon can be de- 
scribed via effective interactions acting on very differ- 
ent length scales: a short-range term describing the ten- 
dency to extend the locally preferred structure and a 
long-range. Coulomb- like term describing the frustration- 
induced free-energy cost associated with this spatial ex- 
tension. Such Coulomb frustrated systems have been 
shown, both by scaling arguments |^ and by Monte 
Carlo simulation^, to display the generic features ob- 
served in fragile glassforming liquids, most notably the 
super-Arrhenius temperature dependence of the relax- 
ation time and the two-step, non-exponential decay of 
the correlation function. 

These Coulomb frustrated models have also been used 
in quite different contexts to describe the formation of 
modulated spatial patterns on mesoscopic length scales, 
such as lamellar and cubic phases in diblock copoly- 
mer meltsB, 0, |6j , microemulsions in water-oil-surfactant 
mixturesMT or stripe phases in high temperature 
superconductors 0. In all these cases, slow relaxation is 
usually observed, and it has been recently argued|jl^ that 
high-temperature superconductors could indeed form a 
"stripe glass" in which glassiness is self generated, i.e., 
does not result from the presence of quenched disorder. 



This latter result has been obtained through an investi- 
gation of the properties of the free-energy landscape of 
the Coulomb frustrated (j)^ scalar field theory: by using a 
thermodynamic approach combining the replica method 
proposed for the study of structural glasses [pd[ [T^ and 
a particular approximation, the self-consistent screening 
approximation (SCSA)|]l^, for calculating the pair corre- 
lation functions, Schmalian and Wolynes [[lO| have derived 
that the free-energy landscape of the Coulomb frustrated 
model becomes non-trivial below some temperature Ta at 
which an exponentially large number of metastable states 
appears; the associated configurational entropy decreases 
with further decrease of the temperature and vanishes at 
a lower temperature T?f [10, 
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Motivated by these results giving evidence for fragile 
glassforming behavior in Coulomb frustrated models, we 
have studied the Langevin dynamics of the Coulomb frus- 
trated 0* scalar field theory within the mode-coupling 
and related approximations. Mode-coupling approaches 
have been widely used to study glassforming liquids |l^, 
pTj , and the dynamical ergodicity-breaking singularity 
predicted to occur in the weakly supercooled liquid re- 
gions, albeit "avoided" in real systems, is taken by many 
as a canonical feature of fragile glassforming systems. 
It is therefore tempting to investigate whether Coulomb 
frustrated models also display this feature. 

The paper is organized as follows. We first present 
the model and summarize the equilibrium phase behav- 
ior and the results previously obtained by Schmalian and 
Wolynes. We also introduce the Langevin equation de- 
scribing the relaxational dynamics of the system. In sec- 
tion III, we derive the evolution equations followed by 
the equilibrium time-dependent correlation function ob- 
tained within two resummation schemes of perturbative 
expansions: the mode-coupling approximation and the 
dynamical SCSA. Section IV is devoted to the search for 
an ergodicity-breaking transition. We find that such a 
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phenomenon is indeed observed with the two approxi- 
mations considered. We also show that the dynamical 
singularity predicted by the dynamical SCSA coincides 
with the temperature Ta at which the replica analy- 
sis of RcfsQ [l^ predicts the occurrence of an expo- 
nentially large number of metastable states. In section 
V, we present the full numerical solution of the mode- 
coupling equations, thereby obtaining the time-evolution 
of the equilibrium correlation function; this latter is sim- 
ilar to that obtained in the idealized mode-coupling the- 
ory of supercooled liquids ||l6| and in mean-field general- 
ized spin-glass models |2^. In the last two section, 



where 0(x) is a real scalar field ((/)k, the associated 
k— Fourier component), V is the volume, u is a strictly 
positive coupling constant, Q is the frustration parame- 
ter, and all momentum integrations are performed up to 
a cut-off A, i.e., |k| < A; ro is a temperature-dependent 
mass which is proportional to the deviation, T — mf^ 
from the mean-field transition temperature of the unfrus- 
trated {Q = 0) model. 

The equilibrium partition function is 

Z = J I?0e-^(^)/^. (3) 

In what follows, we take A = 1, and T is set equal to 1 
in Eq. ^ so that the whole temperature dependence is 
contained in tq . We are interested in the weak-frustration 
region for which Q << 1. 

In the absence of frustration {Q = 0), the model de- 
fined by Eqs. (0-^ reduces to the usual (j)"^ theory. It 
undergoes a second-order transition at a finite temper- 
ature to a broken-symmetry phase characterized by 
a non-zero value of < 0k=o >, where <> denotes an 
equilibrium average. For Q > 0, an ordered phase with 
< ?!'k=o >^ is forbidden, but the system can still un- 
dergo a phase transition at a temperature Tr>o{Q) to a 
phase with long-range modulated order. This transition 
has been studied by Monte Carlo simulation for the case 
of the Coulomb frustrated Ising ferromagnet on a cu- 
bic lattice (2l| and via the self-consistent Hartree approx- 
imation for a Hamiltonian similar to Eq. describing 
microphase separation in diblock copolymer melts it 
has then been shown that, whereas the mean- field the- 
ory predicts a second-order transition, the fluctuations 
change the order of the transition and induce a first- 
order transition. Such a fluctuation-induced first-order 
transition was first discussed by Brazovskii for a related 



we address the question of sensitivity of the results to the 
level and the details of the approximation scheme and we 
give some concluding remarks. 



II. MODEL 

We consider the field-theoretical version of the 3- 
dimensional Coulomb frustrated Ising ferromagnet de- 
fined by the Hamiltonian 

(1) 
(2) 
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FIG. 1: Temperature (ro)-frustration {Q^^"^) equilibrium 
phase diagram in the self-consistent Hartree approximation. 
The full line denotes the fluctuation-induced first-order tran- 
sition to modulated phases; below this line, the paramagnetic 
phase exists in a metastable state. The coupling constant u 
is set equal to 1. 

model Within the self-consistent Hartree approxi- 
mation, the equilibrium (connected) correlation function 

VC{k) 0_k0k > - < 0-k >< 0k > (4) 

is obtained via a self-consistent equation; for instance, in 
the paramagnetic phase where < (/)k >= 0, this equation 
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C(q). (5) 



The renormalized mass, r 
given by 



r — tq + 3m 
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Since fc^ + -p- is minimum for non-zero wave vectors with 
modulus kra = Q^^**, a value characterizing the incipient 
modulated order, one easily checks that r only goes to 
zero when rg — oo. This means that the paramag- 
netic phase is (meta) stable at all finite "temperatures", 
its spinodal being depressed to tq — > —oo. The Hartree 
approximation allows one to calculate the free energy of 
the paramagnetic phase and that of the phase with mod- 



ulated order. One then obtains the temperature 
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of the first-order transition at the point at which the two 
free energies are equal. The details are given in Appendix 
A and the resulting phase diagram is shown in Fig. l|. 

In their recent work, Schmalian and Wolynes |1Q| have 
applied the thermodynamic approach of non-random 
glassforming systems developed by Mezard and Parisi[p^ 
to the Coulomb frustrated model. The basic idea, orig- 
inally motivated by the behavior of a class of mean- 
field generalized spin glasses such as spin and Potts 
glasses is that glassiness arises because of the oc- 

currence of an exponentially large number of metastable 
states. This occurrence, and the associated emergence 
of a non-zero complexity or configurational entropy, 
can be more conveniently studied within the replica 
formalism PI. Approximations are of course necessary 
to solve the corresponding many-body problem and ob- 
tain all relevant correlation functions. Schmalian and 
Wolynes have shown that within the SCSA|i3|, an ap- 
proximation that goes beyond the Hartree result in that 
it explicitly includes more diagrams of the perturbative 
expansion [p3|, there is a temperature Ta at which an ex- 
ponentially large number of metastable states emerges, 
as signaled by a non-zero configurational entropy. The 
configurational entropy decreases with further decay of 
the temperature and it vanishes at a temperature Tk at 
which the system undergoes a random first-order transi- 
tion to an "ideal glass" . 

In the present work, we focus on the dynamics of the 
Coulomb frustrated model defined by Eq. (|l|) . The start- 
ing point is the Langevin equation, 
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dt 
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that describes the purely relaxational dynamics of the 
system; 7]k(i) is a gaussian thermal noise with < 
rii,{t) >= and < ??k(0??k' (i') >= 2T5(k + k')5{t - t'). 
Eq. (0) can be explicitly written as 
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+ ?yk(t). (8) 



Solving this set of coupled non-linear dynamical equa- 
tions is a daunting task, and virtually all available ap- 
proximations amount to performing some self-consistent 
resummation of perturbative expansions, e.g,. expan- 
sions in powers of the coupling constant u or of the in- 
verse of the number of components of the field, 1 /n, for an 
0{n) model. In what follows, we shall consider two such 
self-consistent resummation schemes, the mode-coupling 
approximation and the dynamical SCSA|24|. 



III. DYNAMICAL SELF-CONSISTENT 
APPROXIMATIONS 

To introduce the mode-coupling approximation, we 
first define the time-dependent correlation function 
C(k, t, t') and the associated response function G(k, t, t'): 



(5(k + k')C(k,t,t') = < (/.k(t)0k'(<') >, 
.(k + k')G(k,M')^<|^> 



(9) 



= ^<MtW'(t') > ■ (10) 

As in the preceding section, we set T = 1 in the following. 

The perturbative expansion of C(k, t, t') and (^(k, i, t') 
in powers of u is more conveniently expressed by intro- 
ducing the zeroth-order correlation and response func- 
tions, 
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C7o(k,t,i') ^-^cxp[Ai(k)(f-t')], (11) 

(12) 



>(k) 

Go(k,t,i')=cxp[/i(k)(t-i')] 



where 



M(k) = ro + fc' + ^, 



(13) 



and the two kernels S(k, i,t') and D(k, <,t') defined 
though the standard Dyson equations. 
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C(k,i,t') = f dh f dt2G{k,t,h)[2S{h-t2) + D{k,ti,t2)]G(k,t',t2), 
Jo Jo 

G(k,t,t') = Go(k,<,i')+ / dti /" 'dt2Go(k,i,ti)E(k,ii,t2)G(k,i2,t')- 

Jt' Jt' 



(14) 
(15) 
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The diagrammatic representation of the perturbative 
expansion and a detailed derivation of the modc-couphng 
approximation can be found in Re f . ; the only differ- 
ence with the cases considered in ]19(| is the presence of 
the frustration term Q/k^ in the expression of //(/c), and 
we merely sketch here the main steps of the derivation. 
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The mode-coupling approximation amounts to expand- 
ing the kernels -D(k, t, i') and S(k, i,t') to second or- 
der in u and replacing the bare (zeroth-order) functions 
Go(k, and Go(k, t, i') that appear in the resulting 
expressions by their renormalized counterparts G(k, i, t') 
and G(k,i,t'j. This leads topi: 



DiKt,t')^6u' J ^ J ^G(ki,i,t')G(k2,t,OC(k-ki-k2,i,t'), 

d^ki f (Pk2 
(2^7 (2^ 



(16) 
(17) 



At the same time /i(k) is renormalized to include the 
so-called tadpole diagrams, which replaces Eq. (|l^) by 
an expression similar to that obtained within the static 
Hartree approximation. 



gives 



f dCik,t) 
\ dt 



= -1. 
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In this work, we are interested by the dynamical ^ ^ equation for the response function thus 

properties of the system at equilibrium: therefore, the reads 



fluctuation-dissipation theorem and the time-translation 
invariance apply, which reduces the dependence upon the 
two times t and t' to the mere dependence on the differ- 
ence t — t' and gives 



dGjk, t) 
di 



-/i(k)G(k,i) + / dt'Y,{k,t - t')G{k,t') 
Jo 

(23) 



G(k,i) = -Q{t) 



dC{k, t) 
dt 



and 



s(k,i) = -e(t) 



dD{k, t) 
dt ' 



(18) 



(19) 



with the initial condition G(k, i = 0+) = 1. By Laplace 
transforming Eqs. (|l^), (|l^) and ( |2^ ) and using the ini- 

(?£m\ ^ -G(kA = 0+\ = -I, one 



where 9(t) is the Heaviside step function (8(t) is equal 
to for i < and to 1 pour t > 0). Applying the operator 
Go i(k) = [^(k) + ^] to both sides of Eq. (|l|) yields 



dG{k,t) 
dt 



8{t) - y.(k)G{k,t) 



dt'E(k, t ~ t')G{k, t'), dence leads to 



tial condition (^ "'"^^'^^ j = -G(k,t = 0+) - 
finally obtains 

l = (zG(k,z) + G(k,f = 0)) X 

X (^fi{k) -iz- zb{k, z) - D(k, i = 0)) , (24) 

where G(k, z) = i (iie*^*G(k, t), and a similar expres- 
sion holds for D{k,z). Going back to the time depen- 



(20) 



which, when combined with the time derivative of 
Eq. (@), 



dCjk, t) 
di 



= -{fi{k)^D{k,t^O))G{k,t) 
,,dGik,t') 



I dt'D{k,t-t' 
Jo 



dt' 



(25) 



dG{k,t) 
dt 



e(i) 



t=0 



d^G{k,t) 
dt^ ' 



(21) and (|2 



with the initial condition G{k,t — 0) = 
{fi{k)- D{k,t^O)y^ that follows from Eqs. (|22 
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The mode-coupling approximation finally results in the 
following self-consistent equation for the time-dependent 
correlation function at equilibrium: 



dCQi, t) 

dt 



C(k,i = 0)-iC(k, t) 



[ dt'D(k,t-t') 
Jo 



,,9C(k,t') 



dt' 



(26) 



with 



DiKt) ^6u'J ^ I ^C(ki,t)C(k„t) X 

xC(k-ki-k2,t). (27) 

Except for the absence of inertial term, in 
the purely relaxational dynamics associated with the 
Langevin equation and the cubic dependence of the mem- 
ory kernel D(k, t) on the correlation function C(k, t), the 
above equations are similar to the mode-coupling equa- 
tions used to describe the time-dependent density fluctu- 
ations in supercooled liquids |p6[; they are also analogous 
to those derived for the mean-field spin glass with 4-spin 
interactions |o). 

The necessary input for solving the self-consistent 
Eqs. ( p6| ) and ( p7| ) is the knowledge of the equilibrium 
static correlation function C(k, t = 0). Treating the stat- 
ics and the dynamics of the system on an equal footing, as 
for instance done in the above derivation, leads to consid- 
ering a mode-coupling approximation for the static cor- 
relation function, C{k,t 0) = (At(k) - D{k,t^ 0))"\ 
However, we shall rather introduce more flexibility in 
the mode-coupling scheme (a flexibility that goes with 
the many ways to implement the self-consistency at the 
second order of the perturbative expansion) by allowing 
C(k, t = 0) to be computed with several approximations, 
such as the Hartree approximation and the SCSA, that 
are a priori better behaved than the mode-coupling ap- 
proximation as far as the static properties are concerned. 

A somewhat refined resummation scheme is p rovided 
by the dynamical SCSA. (As mentioned in [ p3| , it con- 
sists in using an n-component vector field cj), resumming 
self-consistently all the diagrams of order 1 /n in the large 
n expansion, and, eventually, for the problem considered 
here, setting n equal to 1.) Details on the derivation of 
the approximate equation for the time-dependent corre- 
lation function C(k, i) can be found in Refs.| 19, |25|| . A 



convenient way to proceed is to introduce a complex aux- 
iliary field cr(x) such that the partition function, Eq. (|^), 
can be rewritten Z = J J V^Vae'"^'^''^^ (here T = 1) 
with 

/ d3x{(V0(x))2+ro^2(^)_^2(^) 



(28) 



2wcr(x)02(x)| 



87r./ ./ x-x' 



The field (7(x) is such that < (t(x) >= ^Ju/2 < (jy^ix) > 
and its connected pair correlation function Co-(x, x') =< 
cr(x)cr(x') >c is equal to 

C,(x,x') =-<5(x-x') 

+u/2{< (j)^{K)(f,^x') > - < (P^x) >< 02(x') >). 

(29) 

One can then apply to the dynamics of the coupled 
fields (t(x) and (/>(x) a treatment similar to that sketched 
above. Defining the equilibrium time-dependent cor- 
relation function Ccr(k, t) via 6(k + k')Ccr(k, i) =< 
o'k(i)ck'(0) >c, the associated kernel D„{'k,t) obtained 
through the Dyson equation (see above), and similar 
functions for the response properties, one can perform 
a mode-coupling approximation to the coupled dynami- 
cal equations for the fields cr(x) and Mx). This leads to 
Eq. (Pq ) with Z?(k, t) now given byPS] 



D(k,t) = 2u 



(2^ 



a(p,OC(k-p,t); 



(30) 



the auxiliary-field correlation function (k, t) is solution 
of the equation 



dt 



-a(k,i-0)-ia(k,t) 

* _,,5a(k,t') 



/ dt'D„{k,t-t')- 
Jo 



dt' 



with 



D„ik,t) = -u J ^C(q,t)C(k-q,i). 



(31) 



(32) 



These equations are supplemented by the initial condi- 
tions C(k, t = 0) = (Ai(k) - L>(k, t = 0)y^ and C<x(k, t = 
0) = — (1 — D„{k, t = 0))^^ which are easily shown to be 
identical to the equilibrium, static SCSA equations first 
derived by Bray| 



IV. TRANSITION FROM ERGODIC TO 
NON-ERGODIC BEHAVIOR 

It is well known that mode-coupling and related ap- 
proximations, when applied to glassforming systems, may 
lead to a dynamical singularity ||l^, ^ . This latter corre- 
sponds to a transition from an ergodic to a non-ergodic 
behavior and is not associated with any thermodynamic 
equilibrium transition. For searching for such a singu- 
larity in the above equations, it is convenient to Laplace 
transform Eq. (p6|), which gives 



C(k, z) 



-C(k,< = 0) 



C{k,t = 0){i + D{k,z)) 



(33) 



An ergodicity-breaking transition is associated with the 
appearance of a non-zero value of the long-time limit of 
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the correlation function, C(k, t — > oo); as a result, in 

the small-z limit, C(k,z) c{k,t-^co) i)(k,z) - 

— -°('^'*~^°°) ^ which when inserted in Eq. (|33|) leads to 



D(k,t 



C(k,t 



C(k, t = 0)(C(k, t = 0)- C(k, t -> c5o)) • 

(34) 



The kernel £)(k, < ^ cx)) is obtained from Eq. (p7|) for 
the mode-coupling approximation and from Eqs. (pOj)- 
for the dynamical SCSA. 



Consider first the equation resulting from the mode- 
coupling approximation, 



C(k,< ^ oo) 



C(k,t = 0)-C(k,i^oo) 



= 6u2C(k,i = 0) 



d^ki f d^k2 



(27r)3 J (27r)3 
I 



C(ki, i oo)C(k2, i oo)C(k - ki - kz, i ^ oo) 



(35) 



Note that this actually represents a set of coupled 
equations for the various k-modes. The necessary input 
for solving this equation is the knowledge of the equilib- 
rium static correlation function C(k,t = 0) (see discus- 
sion above). We consider here two standard approxima- 
tions: 

(i) the Hartree approximation, already presented in 
Eqs. (|) and (|), and 

(ii) the SCSA, described in the previous section and 
leading to 



C (k,0)^Mk) + 2. / (^^^3 i^,n(q) - 



where 



^(k) =ro + k 



2 . Q 



(36) 



(37) 



and 



(38) 



&I, one 



n(k)= J ^C(q,0)C(k-q,0). 

Note that if one neglects the term un(q) in Eq, 
recovers the Hartree approximation. 

In both cases, only the paramagnetic phase < >= 
0, is considered. From now on, we take u = 1 (recall that 
all momentum integrations are cut-off at A = 1). 

Other approximations will be discussed in section VI. 

We have solved the set of coupled equations, Eq. (35), 
by an iterative method, 
two approximations 



We find that with the above 
the mode-coupling approach does 
lead to an ergodicity-breaking transition for the Coulomb 
frustrated model. When decreasing the temperature, i.e., 
the bare mass tq, one reaches a point at which C(k, i — > 
oo) discontinuously jumps to a non-zero value. The tran- 
sition temperature increases as frustration decreases; it 
seems to reach continuously, when Q ^ 0, the (equi- 
librium) critical temperature of the unfrustrated system, 
i.e., of the standard 0"' theory in either the Hartree ap- 
proximation or the SCSA. This behavior is illustrated in 
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FIG. 2: Ergodicity-breaking transition in the ro — (Q^^^) di- 
agram. Dotted line: mode-coupling approximation with the 
static Hartree approximation; dashed line: mode-coupling ap- 
proximation with the static SCSA; full line: dynamical SCSA. 



Fig. ^. One observes a discrepancy between the results 
obtained with the two different approximations (i) and 
(ii) for C(k, t = 0), but it stays within reasonable bounds: 
the relative difference is about 20% or less. As can be 
seen from Fig. ^, the two approximations predict very 
similar correlation functions, both at equilibrium {t — 0) 
and in the non-ergodic state {t — s- oo) when the tem- 
perature is at (or just below) the dynamical transition. 
Roughly speaking, this latter takes place when the maxi- 
mum of the equilibrium correlation function C(k, t = 0), 
a maximum that occurs for |k| ~ km = Q^^^, reaches a 
given, Q-dependent value: this is illustrated for Q — 0.1 
in Fig. ^ where the dynamical transition occurs when 
Maxk{C(k, 0)} ~ 40. As the maximum increases when 
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FIG. 3: Correlation function C(k, t) a.tt = (full line) and in 
the t —* 00 limit (dashed line) in the dynamical mode-coupling 
approximation just below the ergodicity-breaking transition 
for two different values of the frustration Q: (a) static Hartree 
approximation; (b) static SCSA. 



Q decreases slightly more rapidly with the SCSA than 
with the Hartree approximation, the former predicts a 
somewhat higher transition temperature than the latter: 
see Fig. ^. Note that the fact that the ergodicity-breaking 
transition is driven by the maximum of the k-dependent 
equihbrium correlation function is well established in the 
context of mode-coupling approaches jl^ . As illustrated 
in Fig. ^, the transition takes place at the same temper- 
ature for all k-modes. 

The transition form ergodic to non-erogodic dynami- 
cal behavior can also been studied within the dynami- 



FIG. 4: Maximum of the static correlation function 
Maxk{C(k, t)} versus temperature (ro) for t = (full line) 
and t ^ ex (dashed line) in the dynamic mode-coupling ap- 
proximation for a frustration Q = 0.1. The left curves cor- 
respond to the static Hartree approximation and the right 
curves to the static SCSA. The jump in Maxk{C(k, t — > 00)} 
signals the transition from an ergodic (high temperature) 
to a non-ergodic (low temperature) behavior. The inset 
shows the normalized non-ergodicity factor Maxk{C(k, t 
oo)}/Maxk{C(k, t = 0)} versus ro. 



cal SCSA. The corresponding equations to be solved are 
given above and in Appendix A dynamical transition 
is indeed found, and the frustration-dependence of the 
transition temperature is shown in Fig. |2| The predicted 
transition line is not much different from those obtained 
with the above mode-coupling approximations. We show 
in Appendix ^ that the expressions for the correlation 
function C (k, t) for t = and t — > cx) derived within the 
dynamical SCSA when ergodicity is broken are identi- 
cal to those obtained in Refs.[^ |l^ by using the purely 
thermodynamic analysis based on the replica formalism 
and the static SCSA. As a result, the dynamical SCSA 
predicts that the dynamics looses ergodicity precisely 
at the point at which an exponentially large number of 
metastable states occurs in the Schmalian-Wolynes treat- 
ment. Below this temperature, ergodicity is broken: the 
fluctuation-dissipation theorem and the time-translation 
invariance no longer apply, and Eqs. (|2^)-(|3^) should be 
generalized to describe the evolution of two-time corre- 
lation and response functions and the associated aging 
behavior p9|. 

Finally, it is instructive to compare the location of the 
mode-coupling-like dynamical transition with that of the 
equilibrium, thermodynamic transition discussed in sec- 
tion II. The dynamical transition occurs at a tempera- 
ture that is lower than the critical temperature of the un- 
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FIG. 5: Correlation function C(k,t — oo) versus temperature 
(ro) for Q — 0.1 and 3 different momenta k corresponding to 
tfie maximum of tfie function (fcmax, full line), a higher value 
(femax + 0.05, dashed line) and a lower value (fcmax — 0.05, dot- 
ted line). The results are obtained for the dynamical mode- 
coupling approximation supplemented by the static SCSA. 
Note that the jumps occur at the same temperature. 




FIG. 6: Ergodicity-breaking transition (dynamical mode- 
coupling approximation with static Hartree approximation), 
full curve, and fluctuation-induced first-order transition 
(static Hartree approximation), dashed curve, in the rg — Q^^'^ 
diagram. 



V. EVOLUTION WITH TIME OF THE 
CORRELATION FUNCTION 



frustrated system, a temperature that was shown in the 
Monte Carlo study of Ref . ||^ to mark the onset of fragile 
glassforming behavior; it seems to occur at a temperature 
close to that of the fluctuation-induced first-order tran- 
sition from the paramagnetic to the modulated phases. 
This is illustrated in Fig. ^ where we display the first- 
order transition obtained within the Hartree approxima- 
tion (see section II and Fig. |l|) and the dynamical tran- 
sition obtained within the mode-coupling approximation 
supplemented by the static Hartree approximation (as 
discussed above, the other predictions are quite close to 
this latter). 



Actually at small Q's the dynamical transition appears 
even below the temperature of the equilibrium first-order 
transition: in such a case, the dynamical transition takes 
place in the supercooled paramagnetic ("liquid") regime, 
a regime that appears because of the first-order nature of 
the transition to the modulated phases and that can be 
described by the Hartree approximation. As discussed 
in section II, the paramagnetic phase is (meta)-stable 
at all finite temperatures within this approximation, and 
the equilibrium correlation length is therefore finite in all 
the region where the dynamics are studied. 



We have also solved the full set of coupled equations 
describing the time evolution of the equilibrium correla- 
tion function C(k, t) in the mode-coupling approxima- 
tion, Eqs. (|2^)- (|27|). (The algorithm is described in 
Refs[|6[ 1^.) For the input quantity, C{k,t = 0), we 
have used the Hartree approximation. The results are 
shown in Fig. |^ for the time-dependent correlation func- 
tion at a momentum kmax — km = Q^^^ that corresponds 
to the maximum value of the function; curves for the 
frustration parameter Q = 0.1 and several temperatures 
(i.e., several values of the bare mass ro) are shown. One 
observes a behavior typical of the mode-coupling equa- 
tions with a so-called B-type transition[|6j as those used 
to describe glassforming liquids [p^ and those describing 
the dynamics of a class of mean-field generalized spin 
glasses p9[. 

At high temperature, the correlation function decays 
in one step, but as temperature is lowered a second re- 
laxation step appears, that becomes slower and slower so 
that a plateau develops between the two relaxation steps. 
When temperature is further decreased, one reaches a 
point at which the slow ("a") relaxation time diverges. 
The correlation function no longer decays to zero, but 
stays at the plateau value. Below this point, ergodic- 
ity is broken and the mode-coupling equations derived 
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FIG. 7: Time dependence of the correlation function 
C(kmax,t) for Q = 0.1 as predicted by the mode-coupling 
approximation (supplemented by the static Hartree approxi- 
mation). Curves from left to right correspond to decreasing 
temperatures : ro = -1.05, -1.06, -1.065, -1.068, -1.070, 
-1.071, -1.07135,-1.07155,-1.07165, -1.0717, -1.07173, 
-1.071734, -1.0717379, and -1.071738 .The dynamical tran- 
sition is at ro ~ -1.0717379 . . . 



under the condition of equilibrium (with the fluctuation- 
dissipation theorem and the time-translation invariance) 
are no longer valid. 

In the vicinity of the dynamical transition, various scal- 
ing laws are observed, and the slow ( "a" ) relaxation time 
diverges as a power law, 



T^{Q,T)^{T~T,)-\ 



(39) 



where the exponent 7 ~ 1.85 — 1.89 is weakly depen- 
dent on the frustration parameter, provided Q > 0. (For 
(5 = 0, the system shows standard critical slowing down 
with t(T) ~ (T — T'c)^^'', where z is the dynamical expo- 
nent and V the (static) correlation length exponent psf.) 
For illustration, we have plotted the logarithm of the a— 
relaxation time versus temperature rp for two different 
frustrations in Fig. |[ 



VI. SENSITIVITY OF THE RESULTS TO THE 
APPROXIMATION SCHEME 

We have already mentioned (see Ref. ) that enough 
of the non-linearities of the original dynamical equation 
or of the equations in replica space must be kept in any 
approximate treatment in order to find non-trivial phe- 
nomena such as ergodicity breaking and appearance of 
an exponentially large number of metastable states. For 
this reason, the dynamical as well as the replica-space 
Hartree approximations are unable to generate such phe- 
nomena. One must therefore consider improved resum- 
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FIG. 8: Temperature dependence of the a-relaxation time 
Tct obtained from the dynamical mode-coupling prediction 
(with the static Hartree approximation) for the time depen- 
dent correlation function C(kmax, t). Tc is defined as the time 
at which C(kinax,t) ~ 0.1. Two frustrations, Q = 0.1 and 
Q = 0.001, are shown, (a) logj^Q(Ta) versus — ro; the left and 
right curves correspond to Q = 0.1 and Q = 0.001 respec- 
tively, (b) log^Q(rQ) versus logjo(ro - roc(Q)), where roc(Q) 
is the ergodicity-breaking transition point. 



mation schemes such as the mode-coupling approxima- 



tion and the SCSAfM 



The additional point we would like to make here is that 
even in the mode-coupling approximation, the results are 
somewhat dependent upon the supplementary approxi- 
mation which is made to describe the static properties of 
the system entering as an input in the dynamical equa- 
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tion. It is well known, and was recalled above, that the 
location, or even the existence of an ergodicity-breaking 
transition is sensitive to the amplitude of the peak in 
the equilibrium (static) correlation function C(k, t — 0). 
We have seen that the Hartree and SCSA give slightly 
different, but compatible results. If one uses instead a 
somewhat less renormalized version of the SCSA with 
/i(k) in Eq. { ^t\ ) now defined with the Hartree and not 
the full correlation function, i.e.. 



rf3q 1 

(2^Mq) 



(40) 



a different behavior is obtained. As shown in Fig. ^, the 
maximum of the static correlation function appears to 
saturate, as one lowers the temperature, to a value that 
is too small to trigger a breaking of ergodicity. 

Finally, considering the static analog of the mode- 
coupling approximation to compute C(k, t=:0) (see 
Eq. (E3) and below), i.e.. 



C-i(k,t = 0) =ra + k^ + -^ + 3u 
-D(k,t^ 0). 



(2^ 



C(q,i 



= 0) 
(41) 



where D{h, i = 0) is given by Eq. (g^ leads to a situation 
in which the limit of stability (spinodal) of the paramag- 
netic phase is reached at a finite temperature, before the 
occurrence of an ergodicity-breaking transition. 

The validity of these various approximations should of 
course be checked by performing a computer simulation 
of the model. However, one can tentatively conclude from 
the above exercise that despite the formal similarity be- 
tween the dynamics and the statics that comes from using 
the Martin-Siggia-Rose functional formalism! 30 , differ- 



ent levels of approximation may be required to describe 
the dynamical and the static properties of the Coulomb 
frustrated model. 



VII. CONCLUSION 

We have studied the Langevin dynamics of the soft- 
spin, continuum version of the Coulomb frustrated Ising 
ferromagnet. By using the dynamical mode-coupling 
approximation, coupled with reasonable approximations 
for describing the equilibrium static correlation function, 
and the dynamical self-consistent screening approxima- 
tion, we have found that the system's dynamics dis- 
play a transition from ergodic to non-ergodic behavior, 
similar to that obtained in the idealized mode-coupling 
theory of glassforming liquids 1 16 and in the mean-field 
generalized spin-glasses with one-step replica symmetry 
breaking|l^, |l^. This transition occurs in the paramag- 
netic phase, either in the stable or the metastable region. 
It is related to the emergence of an exponentially large 
number of metastable states found by a purely static 
replica approach: the system looses ergodicity because it 
gets trapped in free-energy minima separated from each 




FIG. 9: Maximum of the correlation function C(kmax,t) ver- 
sus temperature (ro) for t = (full lines) and t oo (dashed 
line) in the dynamical mode-coupling approximation. The 
value of the frustration is Q — 0.1. The two upper curves 
correspond to the static Hartree approximation that predicts 
an ergodicity breaking transition. The lower curve that sat- 
urates at low temperatures corresponds to the less renormal- 
ized version of the SCSA (Eq. (0) 
is observed in this case. 



no dynamical transition 



other by infinite barriers. This whole description, as can 
be inferred from the very nature of the approximations 
that amount to partial resummations of perturbative ex- 
pansions, has a mean-field character: thermally activated 
processes are completely ignored. The predicted singu- 
larity is "avoided" in the true dynamics of the system 
and it remains to be seen, e.g., in computer simulations 
of Coulomb frustrated models, what signatures may still 
be observed in the time evolution of the correlation func- 
tion. As for describing the activated processes, other, 
non-perturbative approaches, such as the phenomeno- 
logical frustration-limited domain]^ and entropic-droplet 
pictures |31|, must be used. 



APPENDIX A: FLUCTUATION-INDUCED FIRST 
ORDER TRANSITION 

In this appendix, we calculate the temperature of 
the equilibrium transition between the paramagnetic 
phase and phases with spatially modulated order within 
the self-consistent Hartree approximation. The deriva- 
tion given below closely follows Brazovskii's original 
treatment |^ . 

The starting point is the Hamiltonian in Eq. (^ aug- 
mented by the introduction of spatially varying external 
fields /ik that are linearly coupled to the scalar field 0k- 
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As a result, (/)k is now the sum of an average component, gaussian variational approximation for the fluctuations. 

=< >j and a fluctuation ipk = 4'\s.~ '^k- The self- The resulting equation of state reads 
consistent Hartree approximation is then equivalent to a 



where the connected correlation function C(k, k') ^kV'k' >, is obtained self-consistently by solving 



C-i(k,k')= (ro + fc2 + ^)5(k + k') 



3m 



(27r)3 



(mq77ik+k'-q + C(q, k + k' - q)) , 



r 



(Al) 



(A2) 



together with the inversion formula 



(A3) 



By introducing Eqs. (A4)- (A7) in Eq. ( Al) and recalling 
that km — Q^^"*, one obtains the following equation of 
state: 



In the paramagnetic phase wh en a ll /ik's, and therefore 
TOk's, are equal to zero, Eqs. (A2) and jA^ ) reduce to 
Eq. (|) with C(k,k') = C(k)(5(k + k'). 

In the vicinity of the transition between paramagnetic 
and modulated phases and for small enough frustration Q 
{Q « 1), the modulated order is one-dimensional and 
characterized by a wave- vector km with |km| = km = 
Q^/^. It is then sufficient to consider 



h = (^ro + 

d^k 



+ 3u 



(27r)3r + fc2 



Q 



(r + 2Qi/2 _3„| 



1 m. 



3u|mp 



(A8) 



and 



h^, = h{5{k-k^) + 5{k + k^)) (A4) 



mk = m(<5(k-km) + ,5(k + km)). (A5) 



In this region, the fluctuations of wave-vector k with 
|k| = km are dominant, and Brazovskiip^ has shown 
that the effect of the off-diagonal terms with k 7^ k' could 
be neglected in the correlation function. As a result. 



C(k,k') 



^(k + k!) 

Q 



fc2 



(A6) 



In the Hartree approximation, and below some temper- 
ature, there is a coexistence of the paramagnetic phase 
and the modulated phase. In zero field {h — 0), the for- 
mer is characterized by m = and the latter by m 7^ 0, 
where m is solution of (r + 2Q^I^ — 3w|mp) = 0. The 
transition point, which is then associated with a first- 
order transition, is obtained as the temperature at which 
the free-energies of the two phases are equal. Follow- 
ing Brazovskii[^, it is convenient to calculate directly 
the free-energy difference Ai^(ro) between the modulated 
(to ^ 0) and the paramagnetic (to — 0) phases at a given 
temperature rg from the following expression: 



where the ren ormalized mass in a phase characterized by 
Eqs. (A4,A5), is given by: 



=?'o 




TOkl 



6u|to| 



(A7) 



dF , 
dm'— =2 / dm'h{m') (A9) 
am' 



where hm' is given by Eq. ([A^). One can change the 
integration variable from m' to r' with r' (to') solution of 
Eq. (|A7|). After some algebra, Eq. (A£) can be recast as 



/'r(m) 

uAF= / df 

J r(rh— 0) 



fc2 



Q 
W 



u 
4^ 



dk- 



fc2 



(AlO) 
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where r(m = 0) is solution of Eq. (A7) with m = 0, 
i.e., of Eq. (^), and r( m) a nd ifi are solutions of the two 
coupled equations, Eg. ( A7 ) and (r + 2(3^/^ — 3u|mp) = 
0. By solving Eq. ( AlO ) numerically for several values of 
Q << 1 (and for u = 1), we have found that the sign of 
AF changes at a finite value of ro that marks the first- 
order transition between paramagnetic and modulated 
phases. The result is shown in Fig. |l|. The transition 
being second-order in the mean-field approximation, it is 
then driven first-order by the fluctuations. 



APPENDIX B: ERGODICITY BREAKING IN 
THE DYNAMICAL SCSA 

One can see from Eqs. (^)- ( ^2|) that ergodicity break- 
ing requires that both C(k, t) and (^^(k;, t), and as a con- 
sequence £)(k, t) and -Dcr(k, i], go to non-zero values in 
the limit t ^ oo. From Eq. one obtains 



C(k,i ^ oo) = 



i:>(k,t-^oo)C(k,0)^ 
1 + D{\<i,t-^ oo)C{k,Q)' 



A similar expression can be derived for (k, t 
first Laplace transforming Eq. (|3l|), 



C,(k,z) 



-a(k,o) 



a(k,0)(i-(-Z?,(k,z)) 



(Bl) 
oo) by 

(B2) 



and by looking for the dominant behavior in the small- 
z limit, Ca{k.,z) ~ — Ccr(k, t — > oo)/z, Da-Ck^z) ~ 
—Dcr(k, t oo)/z; one finally gets 



a(k,t 



-jj^(k,t^oo)a(k,o)^ 
i-i?,(k,i^oo)a(k,o)' 



(B3) 



By introducing the time-dependent polarization 



n(k, t)= J ^Ciq, t)C{k - q, t), (B4) 



one can express the memory kernel Dcr(k.,t) given in 
Eq. @ as 



Da(k,t) -■un(k,i). 



(B5) 



so that the t = and i — > oo values of Ccr(k, t) (given be- 
low Eq. ( ]3^ ) and in Eq. (B3), respectively) can be written 
as 



a(k,o) = 

Ca(k, t oo) = 



-1 



1 + un(k, 0) 

-un(k,t ^ oo)c^(k,o)2 
n-un(k,t-> oo)c„(k,o)' 



(B6) 
(B7) 



Recalling that C(k, 0) = (^(k) - D{k,{)))^^ with 
;i(k) = ro + k^ + Q/k^ + u J ^Cp-C(q,0) and that D{k,t) 
zp/en by Eq. (|27|), one obtains with Eqs. (Bl), (B5) and 



a closed set of equations that determines the non- 
ergodicity parameter C(k,t oo). If one changes the 
notations from C(k, 0) and C(k,t oo) to G{k) and 
J^(k), from -D{k,0) and -D{k,t oo) to E(k) and 
Ecr(k), from — uC(k, 0) and —uCa{k,t oo) to I?g(k) 
and I?;(r(k), one can easily check that the above equa- 
tions are identical to those obtained in Refs[p^ Q with 
the replica formalism and the static SCSA. 
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